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Abstract 

The common practice of ignoring the elastic strain gradient in measurements of geomet- 
rically necessary dislocation (GND) density is critically examined. It is concluded that 
the practice may result in substantial errors. Our analysis points to the importance 
of spatial variations of the elastic strain field in relation to its magnitude in inferring 
estimates of dislocation density from measurements. 

1 The role of elastic strain gradients in measurement of 
"geometrically-necessary" dislocation (GND) density 

In the standard operating procedure for measuring (and thinking about) the mesoscale dis- 
location density tensor field, it is generally assumed that the curl of the elastic strain may be 
ignored when the elastic strain field is of small magnitude in a plastically deformed material. 
This belief is widely held among researchers in materials science and mechanics working on 
plasticity and residual stress. For example, a recent published assessment concludes that a 
small lattice strain and its associated gradient field may be neglected in the evaluation of 
the dislocation density, whenever the rotation field is non-uniform. In this short note, we 
demonstrate circumstances in which this conclusion is erroneous. Our discussion depends 
upon the fact that a function of small magnitude may not possess a derivative that is of 
equally small magnitude. 

The main argument is the following. Consider a skew-symmetric tensor field u(xx, x% X3), 
a symmetric tensor field e(xi,X2,x^) and a small nondimensional parameter 0<^<1. We 
assume that the ratio of the magnitude of curl e to that of curl uj is approximately unity or 
less. We define the elastic distortion U e as 
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The Kroner-Nye dislocation density field, denned as 

a = curl U e 

may be approximated by curl uj provided that vcurl e is not of a higher order of magnitude 
than curl uj. For example, let the only non- vanishing components of uj be U12 = —(^21 — kx\. 
Also, let the only non- vanishing component of e be en = kx 2 - Then, the magnitudes 
\curluj\ = k, \curle\ = k, and the only non-vanishing component of ct = curluj + curl {ue) 
is «i3 = k — vk ~ k. However, it is equally clear that if the ratio | curl e\/\curl u>| were to be 
large (e.g. ~ 1 / z^) , then such an approximation does not hold. In what follows, we illustrate 
these remarks with explicit counterexamples. 

2 Micro, and meso-macro scalings 

For simplicity, we consider simply-connected domains for all arguments that follow. Consider 
a symmetric tensor field ey(yi, y 2 , 2/3) as a function of the rectangular Cartesian coordinates 
(yi,i = lto3), where the components of are also expressed with respect to the basis of 
the same coordinate system (as well as all tensor components below). Let us assume that 
the tensor field is twice-differentiable and bounded, and, without loss of generality, that the 
maximum magnitude of any of the components of the field does not exceed unity for all 
values of (yi, y2, 2/3). The field is also assumed to be compatible in the sense of satisfying the 
Saint- Venant compatibility conditions, i.e. it can be expressed as the symmetrized gradient 
(or the strain tensor) of some vector field (the corresponding displacement field): 

d 2 e a _ d 2 e k i _ d 2 e im cF 
dy k dy m dyidy m dy k dyi dyflyi 

We regard this field as an elastic strain field viewed at the microscopic scale. The terms micro, 
meso, macro here refer to the length-scale of variation, and therefore of measurements, of 
various fields. Let us define a small-magnitude, compatible elastic strain field at the 
meso/macroscopic scale according the following definition: 

1 \ (Xi x 2 x z \ 

e ij (x 1 ,x 2 ,x 3 ) := ueij — , — J , (2) 

where < v 1 is the same non-dimensional scale factor introduced previously and m > 2 
is, say, an integer. We use the scaling 

Xi 

v m 

consequently, y^ is a fine (magnified) length scale of observation that resolves the microscopic 
fields of, say, wave-length /, and x^ is a coarse length scale of observation resolving meso- 
macro fields of wave-length L 3> I. We think of v = l/L. 
It is now easy to see that 

d 2 eu d 2 e M d 2 e im | d 2 e km _ 1 / d 2 e u _ d 2 e kl _ d 2 e im d 2 e km \ = Q 
dx k dx m dxidx m dx k dx t dxidxi v 2m - 1 \dy k dy m dyidy m dy k dy t dytdyi J 
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Thus, the elastic strain field is compatible, meaning that there exists a displacement field 
Ui(xi,X2,x 3 ) such that 

1 / dui duj 



2 \dxj dxi 
But then, defining the elastic rotation field as 



Cij. 



1 ( du sl duj 



2 V dxj dxi 



we note that the field, corresponding to the elastic distortion field Uij, vanishes but the 
curl of the elastic rotation field is in general large: 

d 2 Ui dujij <h ij 

U — an — Sikj-z — = — — £ikj~^ r £«£,•■/ t; — 

OXjOXk OXk ox k 

duij f _ N 1 de,ij f xi x 2 x 3 



where are components of the alternating tensor. An example of a displacement field 
Mj whose displacement gradient obviously has vanishing curl but whose strain and rotation 
fields do not is given by U\ = kx 2 x±, u 2 = —kxf/2, u 3 = 0, k a constant. Then, in an 
obvious notation, en = kx 2 and all other strain components are zero, while Ui 2 = kx\ and 
£>13 = 0)23 = 0. 

In the next Section, we explicitly construct strain (e^) and rotation (ttfy-) fields that 
provide counterexamples to the claim that small strains necessarily mean that the strain 
field can be neglected in the determination of the dislocation density field. 



3 Counterexamples 

For the considerations of the previous section to be relevant, it remains to be shown that 
strain fields satisfying (JTJ) in fact exist. But this, of course, is a simple matter. Consider 
any smooth displacement field (2/1? 2/2, 2/3) on a domain, viewed at a resolution of, say, a 
tenth of a micron. By this we mean that when yi are plotted on a scale of a tenth of a 
micron, an unambiguous specification of V{ as a function of these coordinates can be made. 
We require at least that some of the components of the v^s third-order spatial derivatives do 
not vanish. For example, we could take each v^i = lto3, to be a fourth degree polynomial 
in (t/i, 2/2, 2/3) or trigonometric functions. We then define 

[ / Qy ■ Qd ■ \ 

%(2/i,2/2,2/ 3 ) : =2^ + %J (^,2/3). 

By construction, any such field is guaranteed to satisfy (JT]). In terms of this 6^(3/1, 2/2, y 3 ) 
field, define the eij(xi,X2,x 3 ) field from (J2J) for v = 0.01 and m = 2. Roughly, this amounts 
to viewing the same strain field at a scale of a millimeter. Thus, when for example the 
field 6*5(2/1, 2/2, 2/3) contains well-resolved sinusoidal oscillations, the field e^Xx, x 2 , x 3 ) would 
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appear highly oscillatory. The field Uj of the construction above corresponding to this choice 
of 3ij field would be 

3 (x x x 2 x 3 \ 

U i {x 1 ,X 2 , X 3 ) = V Vi ~2 , "a J ■ 

This construction is an explicit example for which t/zere zs no dislocation density in the body 
but if approximated by the elastic rotation gradient, can be erroneously calculated as being 
(arbitrarily) large, depending on the rapidity of the spatial fluctuation of the small-magnitude 
elastic strain field. 

In this counterexample, the dislocation density is exactly zero. However, a small mag- 
nitude, but rapidly oscillating, elastic strain field at the mesoscale would seem to be a 
common occurrence in the presence of a large number of dislocations in the body. Ele- 
ments of the above example show that the curl of the elastic strain in such cases can be 
large in magnitude so that its neglect can lead to substantial underestimation of the ac- 
tual (excess/polar/geometrically- necessary) dislocation density. To illustrate this remark, 
let the elastic distortion be composed of a rotation tensor uj with only non-zero compo- 
nents U12 = — UJ21 = kxi as before, and let the only non-vanishing component of e be 
en = vk sin (x-z/v 2 ), which is a small-magnitude, very high frequency oscillation. Then, the 
only non-vanishing component of ot = curluj + curie is OJ13 = k — (k/u) cos{x 2 / v 2 ), while 
\curluj\ = k. 



4 Concluding remarks 

The discussion so far has been confined to linear kinematics, but even in the geometrically 
nonlinear case we may treat a Right Cauchy Green tensor (RCG) field C = F T F of a smooth 
deformation field with deformation gradient F. Consequently, its Riemann-Christoffel tensor 
is zero, and its essentially unique rotation field R that also satisfies F = R\fC (up to a 
rigid rotation) in general is not curZ-free. Thus, in general, a compatible, i.e. curl-free, F 
field has a non-uniform rotation field with non-vanishing curl. In particular, Shield |Shi73] 
shows that this result holds for compatible finite strain fields of small magnitude with large 
gradients: 

C lk -5 tk = 0(u), C ik ,i = 0(1) 

and 

Rmn.k ^RmpiC pk ,ri ^nk,p) 0(l/). 

Of course, the curl of an RCG field arising from a smooth deformation does not necessarily 
vanish. 

A second remark pertains to the fact that experimental estimates of GND density are 
based on a difference approximation of the derivative where the spatial step-length on the 
mesoscale does not really tend to zero. One may then ask whether the 'discrete' deriva- 
tive, based on a small, but finite, mesoscale step-length Ax, of a continuously differentiable 
function with small magnitude, say, 

9( x ) = u f i£^) , < 1/ < 1, m > 2, 
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is small (we assume / to be a bounded by unity, nondimensional, continuously differentiate 
function). If so, then such a definition would allow the neglect of the small magnitude strain 
field. But even this proposition appears to be untenable. Defining a discrete derivative 
g'(x, Ax) as 

it a \ _ v _ [f iv^ v™) ~ f_ _ v _ \j_ {is™ H™) ~ / \ 

y \ i j a„ Ax 



Ax Ax v m' 

is m 

it can be observed from the middle expression that the term multiplying v is in general not 
going to be small in magnitude since, even for Ax small, the arguments of / are not close. 
Hence, g'(x, Ax) is not going to be small and cannot be ignored. The extreme right-hand-side 
expression shows that for Ax — > 0, < v <C 1 fixed, lim^^o g\ x t Ax) is actually large. 
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